Momentum dependent relaxation rate and pseudogap in doped magnetic insulators 
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The spectral functions and corresponding self energies are calculated within the planar t-t'-J model as rel- 
evant to hole-doped cuprates using the exact diagonalization method at finite temperatures, combined with the 
averaging over twisted boundary conditions. Results show truncated Fermi surface at low doping and t' < 
in the antinodal region while the self energy reveals weakly k- and doping dependent anomalous relaxation 
rate |£"(k, ui)\ ~ a + b\u)\ for u < 0, consistent with recent ARPES results, and a pseudogap-generating 
component of Lorentzian form. The latter is well pronounced at low doping and strongly depends on k and t' . 

PACS numbers: 71.27.-fa, 75.20.-g, 74.72.-h 
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It is well established that high-T c superconductivity in 
cuprates develops directly from a very anomalous normal 
state, which does not follow the usual Fermi-liquid (FL) sce- 
nario and thereby the concept of well defined quasiparti- 
cle (QP) excitations. The nature of the latter is most di- 
rectly accessed by angle resolved photoemission spectroscopy 
(ARPES), which probes the single-electron spectral function 
(SF) A(k, uj) lyj] . Concentrating here on the models for 
hole-doped cuprates, we furtheron discuss several ARPES 
normal-state results, a) Upon doping, the insulating anti- 
ferromagnet (AFM) develops into a metallic paramagnetic 
state where the Fermi surface (FS) appears as a large one, 
consistent approximately with the Luttinger FL volume in 
La 2 _ x Sr 2 ,Cu04 (LSCO) |2] or evidently deviating from it in 
Ca 2 - I Na a; Cu02Cl2 (Na-CCOC) H. b) The FS and corre- 
sponding QP are well pronounced only along the nodal di- 
rection at k ~ (7r/2,7r/2), while FS is truncated or poorly 
defined towards the antinodal points k ~ (ir, 0), where the 
pronounced feature is the large pseudogap first observed in 
optimally-doped Bi 2 Sr2CaCu 2 OK4-,2_(Bi2212) 00] and well 
developed in low-doped LSCO and Na-CCOC H- c) 
The QP relaxation rate, best identified along the nodal direc- 
tion at optimally doped Bi2212 |1, 7], follows the non-FL 
behavior |S"(k, u>)\ ~ b(\co\ + £T), as summarized within 
the marginal FL concept (MFL) |8]. Recently, it has been 
analyzed all along the FS and represented by |S"(k, ~ 
ak + &kM, where cik is strongly momentum dependent and 
large in the antinodal part of the FS, while 6k is nearly a con- 
stant HI. 

Above experimental facts still represent the major chal- 
lenge in the theory of correlated electrons. There are by now 
several numerical results confirming that prototype models 
as the Hubbard model and the t-J model on a square lattice 
can account for such phenomena. Pseudogap behavior within 
the Hubbard model has been observed in the quantum Monte 
Carlo (QMC) studies 1 10], and more recently using dynamical 
cluster approximation (DCA) 1^11^21 and cellular dynamical 
mean-field theory (CDMFT) I13l 1 1 411 revealing its presence 
in both A(k,u) and density of states (DOS) Af(u). Within 
the t-J model as relevant to cuprates, the numerical study 
of SF using the finite temperature Lanczos method (FTLM) 



established the non-FL behavior of S"(k, w) as well as the 
pseudogap in DOS 1151 llfll . More recent studies detected 
also the sensitivity of pseudogap features to the addition of 
next-nearest neighbor hopping t', both in the t-J model fl7Tl 
and the Hubbard model 1 1 lL 1 1 4fl . In spite of the accumulated 
evidence, there is still no consensus on the origin of non-FL 
behavior and pseudogap. However, several numerical studies 
11(1 ll ll Il3i lis. 1 1611 and also analytical approaches 11811 sug- 
gest the short-range AFM fluctuations as the generator. 

In order to distinguish between different scenarios and 
make a comparison with ARPES experiments, more detailed 
knowledge about the momentum dependence of S"(k, ui) is 
needed. In this Letter we present the results for the latter, 
as evaluated within the t-t'-J model, improving the FTLM 
with a continuous k variation on small lattices. We show that 
the self energy can be decomposed for oj < into an essen- 



tially k-independent MFL part |E"(k,w)| 



b\w\ and 



a pseudogap-generating contribution, which is strongly k de- 
pendent. We also find that the pseudogap part essentially de- 
pends on t'. 

Our study is devoted to the extended t-J model 



Sj - -n % n 3 ), (1) 



where cj s = (1 — n^- s )cy is are projected fermionic operators 
not allowing the double occupancy of sites. We include on a 
square lattice besides the nearest- neighbor hopping tij = t 
also the next-nearest-neighbor hopping tij = t' . We consider 
here t' = — Q.3t, as relevant for hole doped cuprates fl7ll . 
as well as the reference t' = 0. In correspondence with ex- 
periments on cuprates we also fix J — 0.3 t and note that 
t - 400 meV. 

We calculate Green's function for projected fermionic op- 
erators, 



G(k,w) = 



dte i ^ t ({c ka (t),4 8 } + ) ) (2) 



and the corresponding SF A(k, oj) = — Im G(k, uj)/tt, where 
/i is the chemical potential. In the following we present re- 
sults obtained by using the exact diagonalization approach for 



small tilted (Pythagorean) square lattices with N = n 2 + m 2 
sites. Although we are primarily interested in the low T — > 
regime, there are several advantages to perform the calcula- 
tion at T > 0, using the FTLM |15]. As discussed before 
fifill . the spectra even on a small system become quite dense 
and macroscopic-like at T > Tf s , where finite-size temper- 
ature depends on the size N as well as on the model used. 
Consequently, the finite-size effects are substantially reduced 
at T > Tf s . For our purpose it is important that in the lat- 
ter regime we are able to extract a meaningful self-energy 
E(k, ui) from the SF. In the following we present results for 
systems of N = 18, 20 sites ranging from low to intermedi- 
ate doping, i.e. Nh — 1,2,3 holes, where the corresponding 
Tf s = 0.1 — 0.15 t. Within FTLM one evaluates separately 
the electron creation and annihilation SF A+ h _(\s., uj) repre- 
senting the transitions N' h — =p 1, respectively. In fact, we 
calculate directly only A^_ h (k, u>) and express A^ h (k, u>) — 
exp(f3uj)A ] ^ h ~ 1 (k,uj). Here we fix /i by requiring the total 
sum rule (1/AOEk/ duA(k,u) =a = (1 + c h )/2. 

On a system with periodic boundary condition (BC) one is 
able to consider SF only for a set of discrete k = ki, I = 1, N. 
In order to scan the whole Brillouin zone, we employ twisted 
BC by introducing the uniform vector potential 9, which mod- 



(31,71)1 (b) 



ifies the hopping elements fjj 



tij = Uj exp(i9 ■ fij) in 



the Hamiltonian, Eq. {I}. It is well known that this proce- 
dure allows one to reach arbitrary momenta k = ki + 9. In 
the following we perform calculation for different phases 9t, 
which are chosen equidistantly in the quarter of the first Bril- 
louin zone (FBZ). In particular, our mesh contains 10 x 10 and 
7x7 k-points for N = 18 and N = 20, respectively. 

Using FTLM we calculate SF and extract corre- 

sponding self energies S(k, oj). Due to non-trivial sum rule 
within the projected model, Eq. Q, the latter are defined via 



G(k,w) 



Ck-s(k, w ) 



(3) 



Since we require that |S(k, u> — » ±oo)| oc 1/uj, the 'free' 
term is uniquely determined as £k = J dujujA(k, uj)/a. It 
should be noted that due to projection the t-J model does not 
directly possess any 'free' band term therefore the dispersion 
Ck is already nontrivial, dependent on doping Ch and model 
parameters 1 1 8 ] . 

In order to reduce finite-size fluctuations of S(k, uj) among 
different k we use additional averaging analogous to the treat- 
ment common in the cluster DMFT approaches. The underly- 
ing idea is that S(k, uj), being quite local quantity in strongly 
correlated systems, varies with k more smoothly than corre- 
sponding A(k, uj). Our results presented furtheron confirm 
this conjecture. Therefore, we perform the averaging of cal- 
culated S(k, uj) (as well as (t) using Gaussian weighting with 
radius 5k ~ 0.3. Such averaged S(k, uj) and £k are then in- 
serted into Eq. l|3} to evaluate SF. 

First, we present in Fig. 1 results for the SF at the chem- 
ical potential, A(k, uj = 0), calculated with t'/t = —0.3 at 
the lowest reachable T = 0.1 1 ~ Tf s , as a continuous scan 




c h =l/20 



Figure 1: SF at the chemical potential A(k, uj — 0) in the first Bril- 
louin zone for the t — t' — J model with t'/t — —0.3. Results are 
shown for T/t — 0.1 and three doping concentrations. Dots show 
points, where SF and E"(k, ui) in Fig. 4 are presented. 



in the quarter of the FBZ. As usual, such a plot is used to lo- 
cate the FS at particular doping. It is well resolved, that at 
lowest doping Ch = 0.05 (Fig. la) we obtain appreciable SF 
weight only in the nodal direction close to k^- = (ir/2,ir/2) 
which is clearly in agreement with experimental results for 
hole-doped cuprates at low doping OlQEIl- Moving to the 
edge of the FBZ one enters the pseudogap region, i.e. a re- 
gion where the QP peaks at the FS are strongly suppressed 
and the FS is in fact not well resolvable, at least not within 
our numerical limitations. While this behavior is most pro- 
nounced at the lowest c/j in Fig. la, it persists also for higher 
but still 'underdoped' c/, = 0.1 in Fig. lb. Evidently, the 
FS gradually builds up with doping, while the pseudogap re- 
mains most pronounced at the edge of the FBZ (antinodal re- 
gion). On the other hand, at 'optimum' doping Ch ~ 0.17 
in Fig. lc the FS becomes almost a continuous line. More- 
over, at the edge of the FBZ near k^ = (tt, 0) the SF be- 
comes even more pronounced than in the nodal direction, the 
effect indeed observed in LSCO and Na-CCOC at intermedi- 
ate doping 1 2, 5]. A decrease in intensity between nodal and 
antinodal region is, however, more reminiscent of electron- 
doped Nd2- x Ce a; Cu04 analogue 1 19]. It should be noted that 
the observed effective FS appears always larger than expected 
from the Luttinger theorem, where Vfs/Vbz = (1 — Ch)/2. 
A deviation is recently established also from the systematic 
ARPES study of Na-CCOC H, whereby our results show 
even larger discrepancy. 
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Figure 2: Same as Fig. 1 for the t — J model. 

For comparison we present in Fig. 2 also corresponding 
plots for the reference t-J model. We detect here more con- 
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tinuous FS line over the entire FBZ already in the low-doping 
regime. Still, the SF weight is much weaker in the latter 
regime, similar to Fig. la. The effect of increasing doping 
is clearly visible in the sharper and better defined FS due to 
the increasing weight of the SF at the FS. The deviation from 
the Luttinger theorem is still observable, although it is smaller 
relative to Fig. 1. The difference between the SF in Figs. 1,2 
can be at lowest doping Ch ~ 0.05 in a naive way explained 
via the rigid band picture and an effective dispersion e(k) of a 
single hole in an AFM 1 20] . While within the t-J model such 
a dispersion is very anisotropic around the k^, leading to a 
minor difference Ae = e(kj<-) — e(kx), t' < induces more 
isotropic dispersion around k^ and therefore more truncated 
FS. On the other hand, at intermediate doping FS already ap- 
pears as a large one with the dispersion according to the renor- 
malized 'free' band where t' directly influences its curvature. 

One of the clearest manifestations of the pseudogap is the 
DOS Af(w) = (2/N)J2k A (Ku), as presented in Fig. 3, 
both for t'/t = -0.3 and t' = 0. We note that the results 
for t-J model are quite close to previous ones 1 16], obtained 
without the BC averaging. This confirms that the pseudogap, 
more or less pronounced at u> ~ 0, is a robust feature of the 
model. Akeady from DOS one can conclude, that pseudogap 
vanishes with increased doping while t' < enhances it. E.g., 
for t' < DOS at u> ~ is reduced due to less coherent band 
and the pseudogap remains more pronounced up to c/, ~ 0.17. 




However, DOS is not selective enough to reveal the origin 
and k-dependence of the pseudogap. Hence, we present in 
Fig. 4 sets of SF and corresponding self-energies S" (k, us) 
taken along the paths shown in Fig. 1 for each doping, re- 
spectively. We choose points, being slightly below the FS in 
the nodal part (as this does not lead to an essential change in 
E"(k, u))) and matching the pseudogapped FS at the edge of 
the FBZ. 

Fig. 4a shows SF A(k,u) and E"(k,w) for t'/t = -0.3 
at lowest doping Ch = 0.05. We notice that along the cho- 
sen path a sharper peak in the nodal region, corresponding 



approximately to a brighter spot of large weight in Fig. 1 a, de- 
velops into a pronounced two peak structure in the antinodal 
region as a manifestation of the pseudogap in SF. Besides that, 
SF shows well known incoherent and nondispersive part for 
u) -c and in this respect strong asymmetry between co < 
and u> > regions, well evident also in Af(u>) in Fig. 3. The 
latter features persist up to high doping. In addition, from 
Fig. 4 it follows that the pseudogap splitting near kx closes 
on doping, whereby the QP become better defined although 
still heavily damped. 




Figure 4: Spectral function A(k, u) (left column) and corresponding 
E" (k, lu) (right column) for t'/t = —0.3 taken in k-points shown in 
Fig. 1 for dopings cj, = 0.05, 0.1, 0.17, respectively. 

A sensitive probe of QP damping and pseudogap behavior 
are self energies E"(k, u>) in Fig. 4 revealing several features: 
a) at uj <C large and broad damping function |E"| 3> t 
is consistent with the incoherent part of SF, b) much smaller 
damping for u> 3> 0, c) linear behavior E" (k, ui < 0) as de- 
scribed within the MFL scenario, d) a pronounced peak of 
a Lorentzian form is building-up at oj ~ as we move to 
the antinodal part of FS, representing a clear signature of the 
pseudogap. Namely, a sharp peak oc itA 2 8(uj) would 
open a real gap in the SF at the FS, while broader one opens 
the pseudogap, with the effective peak splitting Epc ~ 2A. 

It follows from Fig. 4 that the pseudogap contribution in 
S"(k, oj) for t'/t = —0.3 is largest in the antinodal region and 
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Figure 5: Momentum dependence of the pseudogap weight A k ob- 
tained for t'/t = —0.3 and t' = for different along the paths 
presented in Figs. 1,2. 



for lowest doping Ch = 0.05, 0.1, where it is also very sharp. 
Elsewhere it loses weight and becomes indistinguishable from 
the background. In order to quantify S"(k, uj) close to FS and 
to study its momentum dependence we assume for to < the 
form 



E(k,w) ~ T, MFL (k,u) 



At 



/(CO 



ir k ). (4) 



whereby the MFL part is taken as — Y,'l IFL (k, u> ~ 0) ~ ak + 
6k M- We use Eq. (0} to fit the calculated S"(k, u) within the 
range —t < oj < 0.2t. The latter is chosen in this way since 
at u> < the linear part extends at least down to ui ~ — t 
(at c/i = 0.17 even to uj ~ — 2i) while the position of the 
peak oj£ < 0.22; and the linear part at w > is also very 
restricted. Moreover, the fit and the pseudogap interpretation 
is meaningful only as far as I\ < Ak. 

The fitting analysis using Eq. l|4} yields results presented 
in Fig. 4. Again, the variation along the paths shown in 
Figs. 1,2 is considered, whereby only points with Tk < Ak 
are presented. We observe for t' /t = —0.3, as expected from 
Fig. 3, a clear decrease of Ak towards the nodal direction for 
c/i = 0.05, 0.1, whereas for 'optimum' Ch = 0.17 the weight 
is reduced in particular in the antinodal region with a maxi- 
mum at an intermediate k along the path. The t-J model re- 
veals smaller Ak, in particular it is hardly resolved for larger 
Ch — 0.11, 0.17, hence we omit these data in Fig. 4. 

As expected, 6k shows much less pronounced momentum 
and even doping dependence. The linear regime is particularly 
broad for the 'optimum' doping Ch = 0.17, where we find 
6 k = 1-45 ± 0.2, 1.2 ± 0.2 for t'/t = -0.3 and t' = 0, 
respectively. In order to make a quantitative comparison with 
experimentally determined 6 ~ 0.8 0] one can compare only 
the effective damping Zk|S"(k, uj)\ whereby Z^ is the QP 
weight well defined only in the region without pseudogap. In 
ARPES, Zk is less clearly defined and is taken ~ 0.5 |9], 
while we extract Z^ ~ 0.3 ± 0.1. This brings model 6k and 
experimental 6 close together. 

Similar results we get for lower ct < 0.1, i.e. 6k = 
1.8 ± 0.2, 1.5 ± 0.2 for t'/t = -0.3 and t' = 0, respectively. 



However, the range of linearity shrinks for low doping com- 
patible with less justified MFL form for uo ~ 0. Furtheron, the 
MFL scenario predicts a ~ nbT |8]. As our results are calcu- 
lated at T = 0.1 1, we expect the MFL contribution a ~ 0.5 t. 
For the intermediate doping q, = 0.17 and t' = —0.3 1 we in- 
deed find ak close to this value along the nodal direction, with 
the increasing tendency towards the edge of the FBZ where 
ak ~ t. This is consistent with the analysis of ARPES in 
B12212 where in the antinodal region ak ~ 200 meV has been 
observed @|. However, within our approach a more detailed 
analysis of ak is difficult since it is hard to separate it from the 
pseudogap contribution, in particular at lower doping, where 
the pseudogap part dominates. Furthermore, our results in- 
dicate that ak and the pseudogap contribution could become 
the same feature when T > A. Similar question can arise 
also in the interpretation of Bi2212 data where the pseudogap 
appears to be present as well |9]. 

In conclusion, we have shown that within the t-t'-J model 
as relevant for hole-doped cuprates the self energy S(k, u>) 
can be decomposed into two parts. Emfl (k, oj) is almost 
k-independent as well as weakly changing with doping and 
t', but still anomalous and following the MFL scenario. An- 
other contribution generating the pseudogap is strongly k- 
dependent most pronounced in the antinodal part of the FBZ 
and at low doping. It is much more expressed for t' < case, 
which is more relevant for hole-doped cuprates. It seems that 
at intermediate doping the pseudogap contribution transforms 
into a constant part of QP damping ak which therefore be- 
comes k-dependent. 

It should be noted that the pseudogap scale discussed here 
is the large (high-energy) one 1 1 ] as e.g. very evident in LSCO 
J1H and Na-CCOC |J at low doping. The latter is respon- 
sible for effective truncation of the FS in this regime, as ob- 
served in our results in Figs. 1,2. We also emphasize that our 
analysis as well as experimental does not give a clear answer 
whether such a pseudogap indeed vanishes in the nodal direc- 
tion, or the weak coherent QP peak appears inside the pseu- 
dogap. As far as the origin of these phenomena is concerned, 
k and t 1 dependence of the pseudogap as well as underlying 
^mfl(w) are in favor of the interpretation in terms of cou- 
pling to short-range AFM spin fluctuations 
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